Low-energy muon-transfer reaction from hydrogen isotopes to helium isotopes 
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Direct muon transfer in low-energy collisions of the muonic hydrogen H M and helium (He ++ ) is 
considered in a three-body quantum-mechanical framework of coordinate-space integro-differential 
Faddeev-Hahn-type equations within two- and six-state close coupling approximations. The final- 
state Coulomb interaction is treated without any approximation employing appropriate Coulomb 
waves in the final state. The present results agree reasonably well with previous semiclassical 
calculations. 
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Experimental investigation of the low-energy muon 
transfer in collision of muonic hydrogen H M (bound state 
of a hydrogen isotope and muon p~) and He is important 
for muon catalyzed fusion cycle of hydrogen isotopes |Q . 
If the hydrogen is contaminated by He, then the muon- 
catalyzed fusion of hydrogen isotopes will be affected by 
muon transfer from H M to He. This makes the study of 
muon transfer from H^ to He of special relevance. The 
study of such collisions involving three charged particles 
is also interesting from a theoretical point of view as an 
example of rearrangement scattering with Coulomb in- 
teraction in the final state. Such reactions with post- 
collision Coulomb interaction between clusters appear 
frequently in atomic and molecular physics. 

Recently, we formulated S a few-body quantum- 
mechanical description of direct muon-transfer reaction 
from H M to different nuclei using the Faddeev-Hahn-type 
few-body equations |j[ H in the close-coupling approxi- 
mation. We applied this formulation successfully to the 
study of muon transfer among various hydrogen isotopes. 
In that study there was no final-state Coulomb interac- 
tion. Later on we extended this approach to the study 
of muon transfer from hydrogen isotopes to He and Li, 
where the final-state Coulomb interaction is appropri- 
ately included in the calculational scheme H. Next we 
included heavier targets such as, C, O || ||, Ne, S, and Ar 
j7j, in the scheme, where appropriate polarization poten- 
tials were included in addition to the final-state Coulomb 
interaction. In case of these heavy targets the muon 
transfer rates from H„ agree well with the existing ex- 
perimental results H o, Pi. In case of these heavier tar- 
gets there is evidence M pi that the muon is absorbed in 
a specific final state of the heavy nuclei, and hence, the 
muon-transfer rates were calculated in a two-state close- 
coupling approximation including a single state each of 
the target and the incident H^. 

In this report, we perform a two- and six-state calcula- 
tion of direct muon-transfer from H M to He isotopes and 
compare with other calculations. The reaction we study 



present study is meant to test the accuracy and the con- 
vergence trend of our approach. For the lighter target 
such as He, muon transfer could take place via direct as 
well as compound molecular formation. The experimen- 
tal rates include the molecular formation |J . The method 
we consider here is not applicable to muon transfer via 
molecular formation. The experimental transfer rates via 
molecular formation leads to much higher results in the 
case of target He while compared to the theoretical re- 
sult of direct muon transfer. For heavier targets, such 
as C, O, Ne, Ar etc., the molecular rates are highly sup- 
pressed and the muon transfer takes place essentially via 
the direct channel and our theoretical direct rates are in 
agreement with experiment j^, ^, [?[] . 

For the theoretical treatment of a three-body muon- 
transfer rearrangement process, Faddeev-type equations 
H , especially the modified version proposed by Hahn [EJ , 
appear to be very suitable. The two possible asymptotic 
two-cluster configurations of the above rearrangement 
problem are conveniently tackled by a set of two coupled 
Faddeev-Hahn-type equations for components 'Fi and ^ 
of the wave function $ = ^1 + ^2, where each component 
carrys the asymptotic boundary condition for a specific 
configuration |l0| . These equations are very useful to 
incorporate distortion potentials for specific initial and 
final asymptotic states pT| . It is possible to include the 
final-state Coulomb interaction explicitly in these equa- 
tions, so that a low-order approximation to these equa- 



tions produces the correct asymptotic behavior 1 1 1 



(H M )i, + He++ - (He„)+ + H+, (0.1) 

where He ++ represents the He nucleus and (He M )^~ s the 
hydrogen-like bound state of He ++ and a muon. The 



Here we quote the Faddeev-Hahn-type two-component 
equations used in the close coupling type approximation 
for the calculation of muon transfer rates g, § ■ We de- 
note He ++ by 1, the hydrogen isotope (s) by 2 and muon 
by 3. Below the three-body breakup threshold, follow- 
ing two-cluster asymptotic configurations are possible in 
the system 123: (23) — 1 and (13) — 2. These two 
configurations correspond to two distinct physical chan- 
nels, denoted by 1 and 2, respectively. These configura- 
tions are determined by the Jacobi coordinates (7^3, pk): 
= fa - n, p 2 = (f 3 + mifi)/(l + mi ) - r 2 , 
7^3 = F3 - 7^2, pi = (7*3 + mvfty/il + 7712) - fi, where 
fi, rrii (i = 1,2,3,) are coordinates and masses of the 
particle i, respectively. 

Let us introduce the total three-body wave function as 
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a sum of two components Q 

^{ri,r 2 ,r 3 ) = tf 1(7=23, ft ) + * 2 (ri 3 ,p 2 ) (0.2) 

where ^1 (7=23, ft) is quadratically integrable over the 
variable and ^2(^13,^2) over fi 3 . The components 
^1 and carry the asymptotic boundary condition for 
channels 1 and 2, respectively. The second component 
is responsible for pure Coulomb interaction in the final 
state. These components satisfy the following set of two 
coupled equations Q 

[E - (H a + V 23 (f23))]^i(r23,Pi) 

= [(^23(^23) + Vi 2 (ri 2 )) - J7c(p2)]*2(ri3,P2()0.3) 



[E - (H + Vi3(ri 3 )) - U c (P2)}^2(r 13, P2) 

= [(Visits) + ^i2(ri 2 ))]*i(r 2 3, Pi) (0.4) 

where E is the center-of-mass energy, Hq the total ki- 
netic energy operator, Vtj(fij) the pair potential (i ^ 
j = 1,2,3), Uc the final-state Cou lom b in teraction. For 
both Uc = and Uc ^ 0, Eqs. 0.3 and 0.4 are equivalent 



in these cases. The rates A^ r /10 6 sec -1 for both Uc = 
and Uc ^ at low energies are presented in Tabic I 
together with the results of Refs. 13, 14, [l5|. Next 



to the underlying Schrodinger equation for the problem. 
Hence a converged solution of these equations should lead 
to the true dynamical solution of the problem. How- 
ever, an approximate (low-order) close-coupling solution 
of these equations with the proper final-state Coulomb 
interaction Uc should lead to a more converged solution 
when compared with a similar solution with Uc — as 
we shall see in the following. At higher energies the effect 
of the explicit inclusion of the Coulomb potential U c in 
the equations above is reduced and the transfer rates for 
U c = tends towards those for U c ^ 0. 

We solve the integro-differential form of the Faddeev- 
Hahn equation by the close-coupling approximation 
scheme involving up to six states. This procedure consists 
in expanding the wave function components ^1 and ^2 in 
terms of eigenf unctions of subsystem Hamiltonians in ini- 
tial and final channels, respectively. Although, these sub- 
system eigenfunctions are not orthogonal to each other, 
the components ^1 and ^2 satisfy a coupled set of equa- 
tions incorporating the correct asymptotic behavior of 
the wave function. Consequently, there is no problem 
of overcompleteness as encountered in similar expansion 
approaches for rearrangement reactions based on the 
Schrodinger equation. The resultant coupled Faddeev- 
Hahn-type equations are then projected on the expan- 
sion functions. After a partial-wave projection this leads 
to a set of one-dimensional coupled integro-differential 
equations for the expansion coefficients, which is solved 
numerically. The mathematical details of the approach 
have appeared elsewhere and we refer the interested read- 
ers to the original references |2| [| . 

First, we restrict ourselves to a two-level approxima- 
tion by choosing in the relevant close-coupling expan- 
sion the hydrogen-like ground states (H M )i s and (He M )^, 
where H = X H + and 2 H+, and He = 3 He ++ and 4 He ++ . 
Numerically, stable and converged results were obtained 



we extend the calculation to the six-state model where 
we include the H M (ls,2s,2p) and (He p )+(ls,2s,2p) states 
of both muonic hydrogen and helium. These transition 
rates for both Uc — and [/c / are also shown in 
Table I. All rates have converged to the precision shown 
in this table. The present results are consistent with the 
phenomenological isotope effect, e.g., the rate decreases 
from *H to 2 H g2). The effect of including the (2s,2p) 
states in the calculational scheme is also explicit. 

The present six-state calculation, when compared with 
the two-state calculation, shows the trend of convergence. 
We see that at high energies the U c ^ rates tend to- 
wards the U c = rates. This is illustrated in Table I for 
the two-state model for two systems. However, the effect 
of the inclusion of U c is significant at low energies. The 
two-state rates with U c = are the least converged. The 
six-state rates with U c — show the direction of conver- 
gence. From an analysis of our results reported in Table I, 
we find, as expected, the six-state results with U c ^ are 
the most converged. The two-state rates with U c ^ are 
more converged than the three-state rates with U c = 0. 
The explicit inclusion of the final-state Coulomb inter- 
action Uc has enhanced the convergence and has led to 
more consistent rates. The present six-state results for 
Uc 7^ are very close to the semiclassical results of Ref. 
jl3| for all four systems considered. However, the agree- 
ment with the results of Ref. is poor; the present 
rates for 1 H M ( 2 H U ) arc roughly double (ten times) the 
rates of Ref. |L4j . Considering the very qualitative and 
semiclassical nature of the calculation of Ref. 15 , the 
fair agreement with the present few-body quantum treat- 
ment and with Ref. |R| is encouraging. 

The study of three-body charge transfer reactions with 
Coulomb repulsion in the final state has been the sub- 
ject of this report. We have studied such reactions em- 
ploying a detailed few-body description of the rearrange- 
ment scattering problem by solving the Faddeev-Hahn- 
type equations |3, in coordinate space. To provide cor- 
rect asymptotic form in the final state the pure Coulomb 
interaction has been incorporated directly into the equa- 
tions. It is shown that within this formalism, the appli- 
cation of a close-coupling-type ansatz leads to satisfac- 
tory results already in low-order approximations for di- 
rect muon-transfer reactions between hydrogen isotopes 
and light nuclei He ++ . Because of computational diffi- 
culties, in this application we have considered up to six 
states in the expansion scheme (ls,2s,2p on each center 
— Hp and He+), which may not always be adequate. 
Further calculations with larger basis sets are needed to 
obtain accurate converged results. However, the inclu- 
sion of three basis states on each center is expected to 
build in a satisfactory account of the polarization poten- 
tial in the model and should lead to physically acceptable 
results. The agreement of the present rates with those of 
the semiclassical model also assuring. 
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Table I. Low-energy muon-transfer rates A^ r /10 6 sec -1 from 
proton ( 1 H fJ )i s and deuteron ( 2 H M )i s to hydrogen-like 
ground states ( 3 He M ) f g and ( 4 He M ) + s within six-state 
close-coupling model. 



System 


Energy 


Two-state result 


Six-state result 


Other result 




(eV) 


U c = 


U c7 ^0 


U c = o 


u c ^o 




% + 3 He++ 


0.001 


2.1 


8.4 


3.8 


13.6 






0.01 


2.1 


8.4 


3.4 


13.3 


10.9|H.6.3O,7.25lH 

J L ' J II • \i U 




0.04 


2.1 


8.4 


3.3 


12.5 




0.1 


2.0 


8.3 


3.2 


11.0 






5 


1.4 


5.6 










60 


0.5 


1.0 








2 H M + 3 He++ 


0.001 


1.3 


5.2 


2.8 


9.7 






0.01 


1.3 


5.2 


2.7 


9.5 






0.04 


1.3 


5.2 


2.6 


9.4 


9.60, 1.3G3l, 4.770 

l| | i IIJ j i Lr___ll 




0.1 


1.2 


5.1 


2.3 


8.8 




5 


0.6 


3.2 










60 


0.1 


0.15 








1 H M + 4 He++ 


0.001 


1.4 


6.8 


3.2 


12.8 






0.01 


1.4 


6.8 


2.8 


12.4 






0.04 


1.4 


6.8 


2.6 


11.8 


10.711,5.50,6.6511 




0.1 


1.3 


6.7 


2.3 


10.6 


2 H M + 4 He++ 


0.001 


0.8 


5.0 


1.7 


9.6 






0.01 


0.8 


5.0 


1.6 


9.5 






0.04 


0.8 


5.0 


1.6 


9.3 


9.6@,1.0@,4.17|l|] 




0.1 


0.8 


4.9 


1.2 


8.8 





